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Abstract — This paper studies the delay constrained multicast 
capacity of large scale mobile ad hoc networks (MANETs). 
We consider a MANET consists of n s multicast sessions. Each 
multicast session has one source and p destinations. The wireless 
mobiles move according to a two-dimensional i.i.d. mobility 
model. Each source sends identical information to the p des- 
tinations in its multicast session, and the information is required 
to be delivered to all the p destinations within D time-slots. 
Given the delay constraint D, we first prove that the capacity 
per multicast session is O ^min jl, (log p) (log in s p))J^\ j Q 
We then propose a joint coding/scheduling algorithm achieving 
a throughput of O ^min jl, \ ■ O ur simulations show that 

the joint coding/scheduling algorithm achieves a throughput of 
the same order (0 ^min jl, y^})) under random walk model 
and random waypoint model. 

I. Introduction 

Wireless technology has provided an infrastructure-free and 
fast-deployable method to establish communication, and has 
inspired many emerging networks including mobile ad hoc 
networks (MANETs), which has broad potential applications 
in personal area networks, emergency/rescue operations, and 
military battlefield applications. For example, the ZebraNet [1] 
is an MANET used to monitor and study animal migrations 
and inter-species interactions, where each zebra is equipped 
with an wireless antenna and pairwise communication is 
used to transmit data when two zebras are close to each 
other. Another example is the mobile-phone mesh network 
proposed by TerraNet AB (a Swedish company) [2], where 
the participated mobile phones form a mesh network and can 
talk to each other without using the cell infrastructure. 

Despite the importance of these emerging applications, the 
practical deployment of MANETs has been stunned by the 
lack of basic understanding of MANETs. Over the past few 
years, there have been a lot of interest in characterizing the 
capacity of MANETs under a range of mobility models [3], 
[4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], 
[17], [18], [19], [20]. Most of these work assumes unicast traf- 
fic flows and studies the unicast capacity. However, multicast 

'Given non-negative functions /(n) and g(n): f(n) = 0(g(n)) means 
there exist positive constants c and m such that f(n) < cg(n) for all 
n > m; f(n) = Q(g(n)) means there exist positive constants c and m 
such that f(n) > cg(n) for all n > m; f(n) = 0(g(n)) means that 
both f{n) = n(g(n)) and f(n) = 0(g(n)) hold; f(n) = o(g(n)) 
means that linin^co f(n)/g(n) = 0; and f(n) = ui(g(n)) means that 
limn^oo g(n)/f{n) = 0. 



flows are expected to be predominant in many of emerging 
applications. For example, in battlefield networks, commands 
need to be broadcast in the network or sent to a specific group 
of soldiers. In a wireless video conference, the video needs to 
be sent to all the people attending the conference. To support 
these emerging applications, it is imperative to have a fun- 
damental understanding of the multicast capacity of wireless 
networks. In [21], [22]. the authors show that the multicast 

capacity of static ad hoc networks is O ( 1 I per 
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multicast session. In [23], the multicast capacity of delay 
tolerant networks without delay constraints is studied, and 
then the delay associated with the maximum capacity is 
characterized. In [24], the multicast capacity and delay tradeoff 
is established under a specific routing/scheduling algorithm. 
In this paper, we study the multicast capacity of large-scale 
MANETs under a general delay constraint D. We first obtain 
an upper bound on the delay constrained multicast capacity, 
which holds for any communication algorithm. We then pro- 
pose a joint coding/scheduling algorithm with a throughput 
that differs from the upper bound by just a logarithm factor 

In [18], the authors establish the optimal delay constrained 
unicast capacity. The multicast problem differs from the uni- 
cast problem in the following aspects: 

• The capacity of MANETs is highly related to the inter- 
contact rate (the opportunity two mobiles can commu- 
nicate with each other). Since there are multiple des- 
tinations in a multicast session, the inter-contact rates 
between the source and its destinations and the relays 
and their destinations increase. The increase of inter- 
contact rates can improve the capacity of MANETs. On 
the other hand, in the multicast scenario, the information 
needs to be transmitted reliably from the source to all its 
destinations, which generates more traffic in the network 
and requires more transmission resource than that in 
unicast. 

• In MANETs, the mobiles communicate with each other 
using wireless communication. Due to the broadcast 
nature of wireless communication, all mobiles in the 
transmission range of a transmitter can simultaneously 
receive the transmitted packet. In the unicast scenario, 
only the destination of the packet is interested in receiving 
the packet; however, in the multicast scenario, all the 
destinations belonging to the same multicast sessions are 



interested in the packet. Thus, one transmission might 
lead to multiple successful deliveries in the multicast 
scenario, which can increase the capacity of MANETs. 
Due to these differences mentioned, the multicast capacity 
of MANETs obeys a different law from the one for unicast. 

In this paper, we study the delay constrained the delay 
constrained multicast capacity by characterizing the capacity 
scaling law. The scaling approach is introduced in [25], 
and has been intensively used to study the capacity of ad 
hoc networks including both static and mobile networks. We 
consider a MANET consisting of n s multicast sessions. Each 
multicast session has one source and p destinations. The 
wireless mobiles are assumed to move according to a two- 
dimensional independent and identical distributed (2D-i.i.d) 
mobility model. Each source sends identical information to 
the p destinations in its multicast session, and the information 
is required to be delivered to all the p destinations within D 
time-slots. The main contributions of this paper include: 

• Given a delay constraint D, we prove 
that the capacity per multicast session is 
O (min |l, (log p) (log { n sP))\J~^p\^ ■ We then propose 
a joint coding-scheduling algorithm achieving a 
throughput of 8 ^min jl, ypjf^ J ■ The algorithm is 
developed based on an information theoretical approach, 
where a successful delivery can be separated into three 
phases — broadcast, relay and delivery. Each of the 
phase can be modeled as a virtual communication 
channel. Based on the virtual channel representation, 
we propose an algorithm that exploits erasure codes 
to guarantee reliable multicast over the virtual erasure 
channels. The idea of exploiting coding has been used in 
MANETs with unicast flows [16], [15], [18] and mobile 
sensor networks [26]. 

• Finally, we evaluate the performance of our algorithm 
using simulations. We apply the algorithm to the 2D- 
i.i.d. mobility model, random-walk model and random 
waypoint model. The simulations confirm that the results 
obtained form the 2D-i.i.d. model holds for more realistic 
mobility models as well. 

We would like to remark that (a) Similar to the uni- 
cast scenario [3], the mobility significantly improves the 
throughput. While the multicast capacity of a static network 
is O ( ; } ) , our algorithm achieves a throughput of 

\y/n s logn s p J ° ° r 

0(1) when D = n s . (b) Our result again demonstrates the 
substantial benefit of using coding. While the algorithm in 
[24] achieves a throughput of ( — . 1 , ) with an average 
delay Q(\/n s p logp), our algorithm achieves a much higher 
throughput ^^/ pl ° sp ") with the same delay. 

II. Model 

We consider a mobile ad hoc network with n s multicast 
sessions in this paper. Each multicast session consists of one 
source node and p destination nodes as shown in Figure [TJ 
Thus, there are n = n s (p + 1) mobiles in the network. A 
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Fig. 1. A MANET with two multicast sessions, where dst(l,l) and dst(l,2) 
are the destinations of src 1, and dst(2,l) and dst(2,2) are the destinations of 
src 2. A mobile can serve as a relay for other multicast sessions. 

source sends identical information to all its destinations, and 
mobiles not belonging to the multicast session can serve as 
relays. All mobiles are assumed to be positioned in a unit 
torus, where the left and right edges are connected, and top and 
bottom edges are also connected. For the theoretical analysis, 
we assume the mobiles move two-dimensional identical and 
independently distributed mobility model (2D-i.i.d. mobility 
model) [6] such that: (i) at the beginning of each time slot, 
a mobile randomly and uniformly selects a point from the 
unit torus and instantaneously moves to that point; and (ii) 
the positions of mobiles are independent of each other, and 
independent from time slot to time slot. 

Each mobile is equipped with a wireless antenna, and can 
communicate with another mobile within the transmission 
radius. We first assume that each mobile can adapt power 
and use an arbitrary transmission radius, and obtain a general 
upper-bound on the delay-constrained multicast capacity. Then 
we propose a joint coding/scheduling algorithm which (i) 
achieves a near-optimal throughput, and ( ii) requires only two 
transmission ranges {Li,^}, where L\ is for sending out 
information from sources, and is for delivering packets to 
their destinations. 

We also adopt the protocol model introduced in [27] for the 
wireless interference. Let on denote the transmission radius of 
node i, then a transmission from node i to node j is successful 
under the protocol model if and only if the following two 
conditions hold: (i) the distance between nodes i and j is 
less than CKj, and (ii) if mobile k is transmitting at the same 
time, then the distance between node k and node j is at 
least (1 + A)afc (see Figure |2j, where the A > defines 
a guard zone around the transmission. We adopt this protocol 
model because nodes can transmit with different powers (i.e., 
different transmission radius) under this model, which allows 
us to obtain a general upper bound on the multicast capacity 
of MANETs. Note that under this protocol model, the receiver 
of node i associates an exclusion region which is a disk 
with radius and centered at the receiver of node i. 

All exclusion regions associated with successful transmissions 
should be disjoint from each other. We further assume that 
each successful transmission can transmit W bits per time- 



slot. 
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Fig. 4. The virtual channel representation of a multicast session 



Fig. 2. The two transmissions can succeed simultaneously if the distance 
between node j and node k is larger than (1 + AJa^ and the distance between 
node i and node h is larger than (1 + A)a^. 



III. Main Results and Intuition 

In this section, we present the main results of this paper 
along with the key intuition. We use the virtual channel idea 
proposed in [18] to analyze heuristically our system. 
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Fig. 3. The three phases of a typical delivery 

In general, a successful delivery consists of three phases 
(see Figure [3]): 

• Phase-I, the packet is transmitted from the source to some 
relay node; 

• Phase-II, the relay moves to the neighborhood of one of 
the p destinations of the packet; and 

• Phase-Ill, the relay sends the packet to its destination. 
Each of these phases can be thought as a virtual channel as 
in Figure [4] 

« Reliable broadcasting channel: To avoid interference, 
the exclusion regions of successful transmissions should 
be disjoint with each other. To simplify our heuristic 
analysis, we assume all sources use a common transmis- 
sion radius L\ for sending out the information. We also 
assume each exclusion region has an area 7r(ii) 2 .@Here 

2 Note these two assumptions, along with other assumptions introduced in 
this section, are for the purpose of a heuristic argument. Our results hold 
without these assumptions. 



we omit the constant A for simplicity. Thus, the number 
of simultaneous broadcasting at one time slot is at most 
. On average, each source has P\ fraction of time 
to transmit, where 



Pi 



1 



7r(Xi) 2 n s ' 

Thus, the throughput of each broadcasting channel is 

W 

7r(Li) 2 n s ' 

On average, each packet will be received by n(Li) 2 n 
nodes in the neighborhood, and has ir(Li) 2 n duplicate 
copies in the network. 

Unreliable relay channel (erasure channel): We assume 
that all relays use a common transmission radius for 
sending packets to their destinations. The probability that 
a duplicated packets does not fall into the transmission 
range of a specific one of its p destinations during D 
consecutive time slots is 



2\D 



Recall that after sent out from the source, each source 
packet will have 7r(Li) 2 n copies. So the probability that 
none of the duplicated packets falls into the transmission 
ranges of the p destinations during D consecutive time 
slots is 

PndM2= (l-^(L 2 ) 2 )^ (il)2 ", 

which is the erasure probability of the relay channel. 
Reliable receiving channel: Consider the transmissions 
from relays to destinations. When a packet is being trans- 
mitted from a relay, it is delivered to all the destinations 
in the transmission range of the relay. We name one of 
the deliveries as target delivery, and the rest as free-ride 
deliveries. Note that all exclusion regions associated with 
the successful targeted deliveries should be disjoint from 
each other. With a common transmission radius L-2, a 
successful target-delivery associates an exclusion region 
with area ■k(L2) 2 ■ So the number of target deliveries at 
one time slot is no more than 

W 



7T(L 2 ) 2 ' 



Furthermore, along with each target delivery, there are 

( P -1)tt(L 2 ) 2 

free-ride deliveries on average. Thus, we can expect 

W{1 + (p- 1)tt(L 2 ) 2 ) 
tt(L 2 ) 2 

deliveries at each time slot. Since the destinations be- 
longing to the same multicast session request identical 
information, so the throughput per multicast session is 

A + (p - 1)it(L 2 ) 2 W , W(p-1) 



W- 



n s p-K{L 2 ) 2 n s p 



n s pir{L 2 ) 2 
bits per time slot. 
Let A denote the multicast capacity, i.e., the maximum 
throughput per multicast session. Based on the virtual channel 
representation, we can conclude heuristically that 

W 



A = max min 




ir(Li) 2 n s 



where the transmission radii L\ and L 2 solving the maximiza- 



tion are L\ = 9 (-57=) 



and L 2 



e 



We would like to comment that all analysis above is heuris- 
tic, which however captures the key properties determining 
the delay constrained multicast capacity. The rigorous analysis 
will be presented in the rest of the paper, where we will prove 
the following main results: 

Main Result 1: Given the delay constraint D, the multicast 
capacity A (per multicast session) is 



0, if D = 

9(1), ifL» = cj 

o((logp)(\og(n sP )) ' 



(logp) 2 (log(ri a p)) s 



( (lo. 



gp) 2 (log(n s p)) ; 

otherwise . 



Main Result 2: There exists a joint coding/scheduling algo- 
rithm achieving a throughput of 9 (y^J when D is both 

uj(^rT^log(n s p)) and o(n s ). 

IV. Upper Bound 

In this section, we present an upper-bound on the mulitcast 
capacity of MANETs. Note that multicast in MANETs is 
different from unicast in the following aspects: 

• A mobile can send a packet to any of its p destinations, 
which increases inter-contact rates. 

• When a packet is transmitted, it can be received by all the 
destinations in the transmission range, which increases 
the efficiency of the transmission. 

Let Aj[T] to be the number of bits that are delivered to 
destination j before their deadlines expire, up to time T. and 



A P1 = I2j A j[ T }- Furthermore, let B[T] denote the bits 
delivered by target deliveries up to time T. 

Note that in the multicast scenario, one transmission might 
lead to multiple successful deliveries when the destinations 
belonging to the same multicast session are close to each 
other. We first show that the number of occasions that more 
than k(1 + pj 2 )log(n s p) destinations belonging to the same 
sessions are in a disk with radius 7 is small. For a destination 
j, we let H(j,j,t) denote the number of destinations that 
belong to the same multicast session as node j and are within 
a distance of 7 from j at time t. We further define 



(i,7.*)>«;(l+P7 2 ) log(n s p) ■ 



Lemma 1: There exists k > 0, independent of n s and p, 
such that for any 7 € (0, 1] 

T 



E[Z ltK [T}} < 



(1) 



(n s p) 2 
holds. 

Proof: The proof is presented in Appendix A. ■ 
The next lemma establishes a fundamental connection be- 
tween B[T] and A[T]. 

Lemma 2: The following inequality holds: 

E[A[T}] < 5Klog(n s p)E[B[T}] + ^^p(logp) \og(n sP ). 

Proof: The proof is presented in Appendix B. ■ 
Based on the lemma above, we further obtain the following 
two results, which characterize the delay constrained through- 
put without using relays and only using relays, respectively. 

Lemma 3: Consider the 2D-i.i.d. mobility and the protocol 
model. Suppose that packets have to be directly transmitted 
from sources to their destinations. Then, we have 

E[A[T]) < 5 K \og(n s p) I WtJ^^p) 



WkWT 



A 2 



p(logp) log(n s _p). 



(2) 



Proof: The proof is presented in Appendix C. ■ 
Lemma 4: Consider the 2D-i.i.d. mobility. Suppose that 
packets have to be transmitted from relays to their destinations. 
Then, we have 



E[A[T]} < 5K\og(n sP )[ x l—WT(p+l)y/^D 



16kWT 



p(logp) log(n s p). 



(3) 



Proof: The proof is presented in Appendix D. ■ 
The next theorem presents an upper bound on the delay 
constrained multicast capacity. 



Theorem 5: The delay constrained multicast capacity under 
the 2D-i.i.d. mobility and protocol model is 



0, if D = o 

9(1), ifD = «(^ 

o((logp)(\og(n s p)) ' 



(logp) 2 (log(n s p)) 2 
n s 

otherwise 



(4) 



Proof: From Lemma [5] and Lemma |H we can see that the 
throughput by using relay dominates the throughput without 
relay, which implies that the delay constrained multicast 
capacity satisfies: 

E[K[T}} < 5 K \og{n sP )\l — WT [{p+l)^/^D + ^p^ 
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A 2 



p(logp) log(n s p) 



= O 



n s pT (log p)(log(n s p)) J — 
n s 



which leads to the last case 
Note when D = u 



it can be easily 
However, each 



(logp) 2 (log(n s p)) a 

verified that (logp)(log(?7 s p))y^p- = w(l) 
source can send out at most W bits per time-slot, so A < W, 
which leads to the secon,d case 

Next when D = 



(lo S p)"(log(n.p))' J - 11 18 eaS y t0 

verify that DX = o(l). This means under the delay constraint 
D, the information can be transmitted is less than one bit. 
We assume bit is the smallest quantity for information, so the 
capacity is zero in this case. ■ 

V. Joint coding-scheduling algorithm 

In this section, we propose new algorithms that almost 
achieve the upper bound obtained in the previous section. We 
can two different cases: n s = 9(1) and n s = oj(1). For the 
first case, we can use simple round-robin scheduling algorithm 
to achieve the maximum throughput. For the second case, we 
introduce a joint coding-scheduling algorithm that leverages 
erasure-codes and yields a throughput very close to the upper 
bound. 

A. Case 1: n s = 9(1) 

When n s — 9(1), a simple scheme is to let the sources 
broadcast their packets to all the mobiles in the network in a 
round-robin fashion. It is easy to see that the throughput in 
this case is 9(1) per multicast session. 

B. Case 2: n s = lu(1) 

To approach the upper bound obtained in Theorem|5] In this 
subsection, we propose a scheme which exploits coding. This 
scheme achieves a significantly larger throughput than those 
without coding. 

To exploit coding to approach the delay constrained mul- 
ticast capacity, in our algorithm, we code data packets into 
coded packets using rate-less codes — Raptor codes [28]. 



Assume that Q data packets are coded using the Raptor codes. 
The receiver can recover the Q data packets with a high 
probability after it receives any (1 + S)Q distinct coded packets 
[28]. 

We use a modified two-hop algorithm introduced in [3], 
which consists two major phases — broadcasting and receiv- 
ing. At the broadcasting phase, we partition the unit torus into 
square cells (broadcasting cells) with each side of length equal 
to 1/ yfnli which is of the same order as the optimal L* . All 
sources use a transmission radius y/2/ ' ^fn~ s in the broadcast- 
ing phase. To avoid interference caused by transmissions in 
neighboring cells, the cells are scheduled according to the cell 
scheduling algorithm introduced in [25] so that each cell can 
transmit for a constant fraction of time during each time slot, 
and concurrent transmissions do not cause interference. We 
assume each cell can support a transmission of two packets 
during each time slot. In the receiving step, the unit square 
is divided into square cells (receiving cells) with each side of 
length equal to 1/ \Jn s p 2 D. The transmission radius used in 
this phase is \/2/ {Jp 2 n s D. 

Similar as in [18], we define four classes of packets in 
the network: We also define and categorize packets into four 
different types. 

• Data packets: uncoded data packets. 

• Coded packets: Packets generated by Raptor codes. 

• Duplicate packets: Each coded packet could be broadcast 
to other nodes to generate multiple copies, called dupli- 
cate packets. 

• Deliverable packets: Duplicate packets that are in the 
same destination with one of its destinations. 

Joint Coding-Scheduling Algorithm: We group every 2D 
time slots into a supertime slot. At each supertime slot, the 
nodes transmit packets as follows: 

(1) Raptor encoding: Each source takes -^^/Djn s data 
packets, and uses Raptor codes to generate D coded 
packets. 

(2) Broadcasting: This step consists of D time slots. At 
each time slot, in each cell, one source is randomly 
selected to broadcast a coded packet to 9(p + 1)/10 
mobiles in the cell (the packet is sent to all mobiles 
in the cell if the number of mobiles in the cell is less 
than 9(p+ 1)/10). 

(3) Deletion: After the broadcasting phase, all nodes check 
the duplicate packets they have. If more than one du- 
plicate packet belongs to the same multicast session, 
randomly keep one and drop the others. 

(4) Receiving: This step requires D time slots. At each time 
slot, if a cell contains no more than two deliverable 
packets, the deliverable packets are broadcast in the cell; 
otherwise, no node in the cell attempts to transmit. At the 
end of this step, all undelivered packets are dropped. The 
destinations decode the received coded packets using 
Raptor decoding. 

Theorem 6: Suppose D is both uj(^/n^ log(n s p)) and 
o(n s ), and the delay constraint is D. For sufficiently large n s , 



at the end of each super time slot, every source successfully 
transmits 

D ID 
500 1 

packets to all p destinations with a probability 1 — 

Proof: We follow the analysis in [18] to prove the 
following three steps: 

(1) Step 1: During the broadcasting step, with a high 
probability, a source sends out -j coded packets; 

(2) Step 2: After the deletion step, with a high probability, 
a source has at least ^ coded packets that each of them 
has more than ^ duplicate copies in the network. 

(3) Step 3: Each destination receives more than -^y 
distinct coded packets after the broadcasting, which 
guarantees that it can decode the original ^ ^ / -P- data 
packets with a high probability. 



Analysis of step 1 

Let Bi[t] denote the event that node i broadcasts a coded 
packet to 9(p + 1)/10 mobiles at time slot t. According to the 
definition of Bi[t], we have that 

Pr (Bi[i\) = Pr (> 9p/10 mobiles in the cell) 

• Pr (i is selected| > 9p/10 mobiles in the cell 
> Pr (> 9p/10 destinations in the cell) 

• Pr (i is the only source in the cell) . 

Since the nodes are uniformly and randomly positioned, 
from the Chernoff bound, we have 

Pr (> 9p/10 destinations in the cell) > 1 — 2e~ 5°° . 

Note that there are n s sources in the network, so 



of distinct mobiles receiving duplicate packets from session 
i. Assume that source i sends out Df, coded packets. The 
number of duplicate copies left after the deletion is the same 
as the number of nonempty bins of the following balls-and- 
bins problem: There are n s p— 1 bins. At each time slot, 9p/10 
bins are selected to receive a ball in each of them. This process 
is repeated by Db tunes. 

Let Ni to be the number of duplicate packets belonging to 
multicast session i after the deletion. From Lemma 22 in [18], 
we have 

Pr (JVi > (1 - S)(n s p - l)pi) > 1 - 2er s2 ^ p ~ 1 ^ 1 / 3 , 



where 

(n s p - l)pi 



1 — 2e aoo 



1-1 

n. 



which implies that for large p and n s , we have 

Pr(Bi[t]) > 0.36. 

Then from the Chernoff bound again, we have that for suffi- 
ciently large D, 



I 1b<[*] > y ) > 1 - e 3 °°° 



(5) 



Thus, with high probability, more than D/3 coded packets 
are broadcast, and each broadcast generates 9p/10 copies. 

Analysis of step 2 

For analysis purpose, we dropped some of the duplicate 
packets to guarantee that a mobile carries at most one packet 
for each multicast session other than the session it belongs to. 
We next study the number of coded packets that have more 
than 4p/5 duplicate copies. 

Note that the number of duplicate packets of session i 
left in the network after the deletion is equal to the number 
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where the last inequality holds for sufficiently large n s (recall 
that D = o(n s ) under the assumption of the theorem). 
Choosing 5 = 1/50, we have that for sufficiently large n s 
and p, 



Pr 



22 



D 
t=l 



D h > 1 - 2e" 



(6) 



Given that there are more than 22 Di,p/ 25 duplicate packets 
left in the network, we can easily verify that more than 
2£>b/5 coded packets will have Ap/5 duplicate copies because 
otherwise less than 22£>tp/25 duplicate packets would be left. 
Letting A; denote the number of coded packets of session i, 
which has more than 4p/5 duplicate packets after the deletion, 
we have 

D 




D \ „ D_ 

> — > 1 — 2e 10000 . 



(7) 



Note that after the deletion, all duplicate packets belonging 
to the same multicast session are carried by different mobile 
nodes. 

Analysis of step 3 

We consider a coded packet of multicast session i, which 
has at least -j? duplicate copies after the deletion. Let T>i[t] 
denote the event that the coded packet is delivered to its I th 
destination at time slot t. 

First we consider the probability that one of the duplicate 
copies of the coded packet is in the same cell with its I th 
destination. In the receiving phase, we use the cell with each 
side of length equal to 1/ \J n s p 2 D, so the average number of 
nodes in each cell is 



n s (p+ 1) 
\Jn s p 2 D 



Recall that the duplicate packets belonging to the same mul- 
ticast session are carried by distinct mobiles after the deletion, 
so their mobilities are independent. Assuming the number of 
duplicate copies of the coded packet under consideration is 
M, we have 

Pr (only one copy is deliverable to the Z th destination) 

M-l 



M- 



yj n s p 2 D 

Note that M < p, so as n. 



1 - 



\/ n s p 2 D 

oo, we have 
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M-l 



yJn s p 2 D J 
For sufficiently large n s , we have 

Pr (only one copy is deliverable to the I th destination) 

39 



> 



Next, we consider the probability that the duplicate copy is 
delivered given that it is the only copy which is deliverable 
to the I th destination. Suppose we have M nodes in the 
cell containing the I th destination. According to the Chernoff 
bound, we have 




> 1 



(9) 



Note the deliverable copy to the I th destination will be 
delivered if the M — 2 other mobiles (other than the mobile 
carrying the copy and the I th destination for the copy) do not 
carry deliverable packets and there are no deliverable packets 
for the M — 2 mobiles. 

Now given K mobiles already in the cell, we study the 
probability that no more deliverable packet appears when we 
add another mobile. First, the new mobile should not be the 
destination of any duplicate packets already in the cell. Each 
mobile carries at most D duplicate packets, so at most KD 
duplicate packets are already in the cell. Each duplicate packet 
has p destinations. For each duplicate packet, we have 



Pr (the new mobile is its destination) 



Thus, from the union bound, we have 



V 



n s {p+l) - K 



Pr (the new mobile is a new destination) 



< 



pKD 
n s (p+l)-K- 



(10) 



Note that each source sends out no more than D duplicate 
packets and each duplicate packet has no more than p copies, 
so at most KDp mobiles carry the duplicate packets towards 
the K existing mobiles in the cell, and 



Pr (new added mobile brings new deliverable packets) 
KDp 



< 



n s (p+ 1) 



(ID 



From inequalities ( [Tol l and ( TTTb . we can conclude that the 
probability that the new added mobile does not change the 
number of deliverable packets in the cell is greater than 

IKDp 



1 - 



n s (p+l) - K' 



Starting from the mobile carrying the duplicate packet and the 
I th destination of the packet, the probability that the number 
of deliverable packets does not change after adding additional 
M — 2 mobiles is greater than 



M 



n i 



K=2 



2KDp 



n s (p+ 1) - K 



> 1 



2MDp 



n s (p+l)-M 



M-2 



When M < j^^pj^, we have that for sufficiently large n s , 

2MDp 



n s (p+l) - M 



M-2 < 2.5, 



(8) and 



M 

n 

if =2 



1 



2pKD 



n s (p+l)-K 



> 



-2.5 



(12) 



Now according to inequalities ©, ©, and (fT2l . we can 
conclude that for sufficiently large n s , 

1 1 



Pr (©,[*])> 



16 v^Td' 



(13) 



which implies at each time slot, a coded packet with at least 
4p/5 duplicate copies is delivered to its I th destination with a 
probability at least —4 „ . 

1 J 16Vn s D 

Note at each time slot, one destination can receive at most 
one packet. So the number of distinct coded packets delivered 
to the 2 th destination of multicast session i is the same as the 
number of nonempty bins of following balls-and-bins problem: 

Suppose we have ^ bins and one trash can. At each time slot, 
we drop a ball. Each bin receives the ball with probability 
16 V» d 1 ana * tras h can receives the ball with probability 
1 — P, where 

D 



120V"7D' 

Repeat this D times, i.e., D balls are dropped. Note the 
bins represent the distinct coded packets, the balls represent 
successful deliveries, and a ball is dropped in a specific bin 
means the corresponding coded packet is delivered to the 
destination. 

Let Xi i denote the number of distinct coded packets 
delivered to destination I of session i. Under the condition 
that at least 2D /15 coded packets of session i have more than 
4p/5 duplicate copies each, Xij is the same as the number of 
nonempty bins of the above balls-and-bins problem. Choose 
5 = 1/6. From Lemma 22 in [18] we have 



Pr Xu > 



L — 15 



> 1 - 2e 



__D_ / J - 16V" 3 i> 
810 1 1 L 



Using the fact that 1 — e x > x — x 2 /2 for any x > 



Pr Xu > — J — 



2D \ el- /x 

Ai> > 1 - 2e 13000 V »• . (14) 

15 / 



Note that D^J ^- — ► oo under the assumption of the theorem 

(D = oj^n^log(n s p)). 

Summary 

Combining inequalities (0, (0 and JT4t . we can conclude 
that 




__D_ D -D / g 

> 1 — g 3000 — £ 10000 — 2e 13000 V n s 

Furthermore, for sufficiently large n s and p, we also have 
Pr ( X ^^/? fora11 ^) 

/ __D_ D 2_,/~D 

> 1 — n s p e 3000 — e 10000 — 2e 13000 V 

> i--L. 

n s p 

where the last inequality holds under the assumption of the 
theorem (D = u(^/n^log(n s p))). Note that a destination can 

decode the x — data packets after getting * / — coded 

500 y n s * e> & 400 y n s 

packets with a nigh probability, so the theorem holds. ■ 
From the theorem above, we can see that the throughput 
per multicast session is 




VI. Simulations 

In this section, we use simulations to verify our theoretical 
results. We implement the joint coding-scheduling algorithm 
for different mobility models, including 2D-i.i.d. mobility, 
random walk model and random waypoint model. We consider 
an MANET consisting of n s multicast sessions, and the 
mobiles are deployed in an unit square with n s sub-squares. 
The random walk model and random waypoint model are 
defined in the following: 

. Random Walk Model: At the beginning of each time 
slot, a mobile moves from its current sub-square cell 
to one of its eight neighboring sub-squares or stays at 
the current sub-square. Each of the actions occurs with 
probability 1/9. 
. Random Waypoint Model [29]: At the beginning of 
each time slot, a mobile generates a two-dimensional 
vector V = [V x , V y ] , where the values of V x and V y 



are uniformly selected from [1 / ^Jnl, . The mobile 

moves a distance of V x along the horizontal direction, and 
a distance of V y along the vertical direction. 



A. Multicast throughput with different numbers of sessions 

In this simulation, the number of multicast sessions (n s ) 
varies from 200 to 1000, each multicast session contains p 
1 destinations, and the delay constraint is set to be 2D = 200 
time slots. Figure [5] shows the throughput per 2D time slots 
of the three mobility models with different values of n s 

Our theoretical analysis indicates that the throughput is 
9 (2DJ^p) . To verify this , we plot a (2D^^j in 
Figure [5] where a = 0.09 is obtained by using Matlab to fit 
the simulation data of the random walk model. Our simulation 
result shows that the throughput under the three mobility 
models all evolves as 6 (jjljp) ■ Also, the 2D-i.i.d. mobility 
has the largest throughput and the random walk model has 
the smallest throughput. This is because the distance a mobile 
can move within a time slot is the largest under the 2D-i.i.d. 
model and is the smallest under the random walk model. Our 
results indicates that the throughput is an increasing function 
of the mobility speed ( the distance a mobile can move within 
a time slot). 
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Fig. 5. Throughput per multicast session per 2D time slots with different 



B. Multicast throughput with different delay constraints 

In this simulation, we fix n s — 500 and p = 10, and 
change D from 100 to 400 with a step size of 50. We also 
use Matlab to fit the data under the random walk model to get 
the coefficient a = 0.075. For all three mobility models, the 
simulation results match the theoretical order result. 

C. Multicast throughput with different session sizes 

In this simulation, n s = 500, the delay constraint is set to 
be 2D = 200, and p varies from 4 to 40 with a step size of 
4. Figure |7] shows that the throughput is almost invariant with 
respect to p. 

From the simulations above, we can see that the 
throughput is achievable not only under 2D-i.i.d. model, but 

3 In our simulations, we only count the number of distinct packets delivered 
that are successfully delivered before their deadlines expire. We do not 
consider coding and decoding in our simulations. 
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Fig. 6. Throughput per multicast session per 2D time slots with different 
delay constraints 



S 25 



Q 20 



■a 15 



10 



2 

i— 0. 



-e— 2D-i.i.d. 

Random walk 
■ $ - Random waypoint 



° 5 — § — ■© — e — -e — o o — e— 



10 



20 
P 



30 



40 



Fig. 7. 

p's 



Throughput per multicast session per 2D time slots with different 



also under more realistic models such as random walk model 
and random waypoint model as well, which indicates that the 
theoretical results we obtain based on the 2D-i.i.d. mobility 
model hold for more realistic models as well. 

VII. Conclusion 

In this paper, we studied the delay constrained multi- 
cast capacity of large-scale MANETs. We first proved that 
the upper-bound on throughput per multicast session is 

O (min |l, (logp)(log i n sP)) ^/^j) , and then proposed a 
joint coding-scheduling algorithm that achieves a throughput 
of 8 fmin |l, ^/n }) ' ^ e a ^ so va lidated our theoretical 
results using simulations, which indicated that the results based 
on 2D-i.i.d. model are also valid for random walk model and 
random way point model. In our future research, we will 
study ( i) the impact of mobile velocity on the communication 
delay and multicast throughput; and ( ii) the delay constrained 
multicast capacity of MANETs with heterogeneous multicast 
sessions, e.g., different multicast sessions have different sizes 
and different delay constraints. 
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Appendix A: Proof of LemmaQ] 

Recall that each multicast session contains p destinations. 
The probability that a mobile is within a distance of 7 from 
node j is ttj 2 . Thus, H (j, 7, t) is a binomial random variable 
with p — 1 trials and probability of a success irj 2 , and 

E[H{j, 1 ,t)] = {p-l)^ 2 . 

Now choose n such that 

k(1 + p7 2 ) log(n s p) > (p - 1)7T7 2 . 



10 



(15) 



Note that M / p ilY 1 , 7 < n for n s p > 3, so we can choose 

(1+P7 ) log(n s p) b r ^ 1 

k independent of n s and p. Next, define 

= n(l + jry 2 )\og(n s p) _ ^ 

(p — l)ivy 2 

which is positive due to inequality (1151 . 

According to the Chernoff bound [30], we have 

Pr (H(j, 7, t) > k(1 + p 7 2 ) \og(n sP )) 

/ <5 \ (P- 1 ) 7r 7 2 

< 



( P -1)ttj 2 {1+S) 



< 



1 + 5 



1 + 5 



k(1+p7 2 ) log(n s p) 



k(1+P7 ) log(n s p) 



k(1+P7 2 ) log("sP) 
(p-l)7T7 a 

<, , p-re(l+P7 2 ) log(n a p) 
— (a) e 

< g-Klog(" s p) 

where inequality (a) holds for any n such that 



< e 



(16) 



K(l+p7 2 )(logp+logn 3 ) 
(P-1)T7 2 



Thus, we can conclude that there exists n > 0, which is 
independent of n s and p, such that 



E[Z^ K [T}} < E 



^2 ^2 1 H(j, 1 ,t)>K 1 *\og(n s 
j:j is a destination t — 1 



- ^2^2 E (i,7,*)>«7 2 log(n s p)] 

and the theorem holds by guaranteeing k > 2. 

Appendix B: Proof of Lemma[2] 

First, we present some important inequalities that will be 
used to obtain the upper-bound on throughput. Let 1Z[T] 
denote the number of bits that are carried by the mobiles 
other than their sources at time T (including those whose 
deadlines have expired), and as the transmission radius used 
to deliver bit B. The following lemma is presented in [18]. 
Inequality ( fTTI i holds since the total number of bits transmitted 
or received in T time slots cannot exceed n s pWT. Inequality 
( TT8l ) holds since the total number of bits transmitted to relay 
nodes cannot exceed n s (p+l)WT. Inequality $1% holds since 
each successful target delivery associates an exclusion region 
which is a disk with radius Ao.b/2- 

Lemma 7: Under the simplified protocol model, the follow- 
ing inequalities hold: 



API 
\K[T}\ 

B[T] A 2 



< 
< 

< 



n s pWT 
n s {p+l)WT 

WT 



(17) 
(18) 

(19) 



where 



B=l 

\TZ[T}\ 



is 



the cardinality of set 



TZ[T}. 

□ 

We index the target deliveries using B. Let (3b denote the 
number of deliveries associated with target delivery B. Given 
a 7 G [0, 1], we classify the target-deliveries according to q^. 
We say a target-delivery belonging to class (7, to) if 2 r ' 
clb < 2™7. Thus, A[T] can be written as 

~B[T] 

E[h[T}} < E 



1 



< 



<7 



B = l 



\- log 2 7l 



B[T] 

/3Bl2 m - 1 7<a B <2 m 7 

5 = 1 



Note that (3b > k(1 + 4p7 2 ) \og(n s p) implies that 

H(d B , 27, t) > k(1 + 4p 7 2 ) log(n s p) 

as shown in Figure [8] where (is is the destination receiving 
the target delivery, so we have 

(1 + %>7 2 )log(?vp)l a B <7 

H(d B ,2j,t)<K{l+4p- 1 2 )log(n B p 

(20) 



+ (3b 1 a B <7 

-H"(d B ,27,t)>K(l+4p7 2 ) log(ra 3 p) 



# T node./ (type i) 



\ area 1' 



area 2 ' 

Fig. 8. H{d B ,2~f,t) > p B 
Furthermore, it can be easily verified that 

B[T] 

a B <f 



i, and D(si,t) the distance between source Sj and its nearest 
destination, i.e., 

D(si,t) = min dist(si, di j)(t). 
i<j<p 

Thus, we have 

Pr(D( Sl ,t)<L)<l-(l-nL 2 )P, 
which implies 



E 



T n s 



EE 1 



D( Si ,t)<L 



< Tn s pirL 



< 



B=1 ff(d B ,27,t)> K (l+4p 7 -)log(n i ,p) 

r- io s 2 ti 

E * 

m— 1 

pWT 

(a) o o~ 

v ; n%p z 
WT 

9 ' 

7l|p 



Since at most bits a source can send during each transmis- 
sion, we further have 



E [B[T]] = E 



~B[T] 




~B[T] 


E 1<*B<L 


+ E 


E ^a B >L 


B=l 




B=l 



B[T] 

E 2 m - 1 7<Q fl <2 m 7 
B=l ff(d B ,2 m + 1 7,*)>K(l+P7 2 2 2m+2 )log(n s p) 



(21) 
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< WE 






E lo B> L 








5 = 1 



< WTn s pirL + E 



B[T] 

E ^ a B>L 
B = l 



Next, applying the Cauchy-Schwarz inequality to inequality 

where inequality (a) yields from inequality £[6}. ®' we can obtain that 

Now from the inequalities (|20T > and d2"Tl i. we have for any / b[t] \ 2 /b[t] \ /b[t] 

< 7 < 1, 

£7[A[T]] (22) 



B[T] 

E 

,B=1 



as 



< Ei £(«*) s 



B=l 
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n 2 p 



+ Klog(n s p) (l + 4p7 2 )E 



B[T] 

y ] las<7 
B=l 



< BIT 



AWT 

^A 2 "' 



which implies that 



< 



(b) 



[-log 2 7l S[T] 

X; (i+p2 2 ™+v)i? E 1 ^- 

m=l B=l 

K log(n s p)(l + 4p 7 2 )£;[B[r]] 



7<ctB<2 m 7 



ttA 2 



V^BfT]] > E 
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I- log 7/ log 2] 

k log{n s p) ^2 P 2 
WT 



AWT 
ttA 2 

B[T] 

E aB 

B=l 



!m+2^2_ 



WT 



2 2m - 2 A 2 7 2 



+ «logKp)(l + 4 P7 2 )£;[S[T]] 



B[T] 

E ^ a B>L 
B=l 



> LE 

B=l 

> L (-^[S[T]] - WTrisp-nl?) , 



ttA 2 log 2 



(-log7)plog(n s _p) 



< 5k log(ri s p)S [B[T]] + 16 ^f T p(logp) log(n s p), 

where inequality (b) yields from inequality ( fT9] l, and the last 
inequality holds when 7 = -4=. 

Appendix C: Proof of Lemma[3] 

We first bound the total number of targeted deliveries under 
the constraint that sources need to directly send information 
to their destinations. Let s; denote the source of multicast 
session i, d^j denote the j th destination of multicast session 



where the first inequality follows from the Jensen's inequality. 
Now we choose L = \ tM§¥-^ , we can obtain that 

y 2WTn s p7T ' 

WT^^p>E[B[T]]. 
By substituting into the bound on A[T] in Lemma [2 we have 



E[A[T}} < 5 K \og{n s p)[WT x /^^p- 



16kWT 



A 2 



p(logp) log(n s p). 



Appendix D: Proof of LemmaO 

Denote by H(b) the minimum distance between the relay 
carrying bit b and any of the p destinations of the bit during 
D consecutive time slots. We have 

Pr(if(6) <L) < l-(l-nL 2 ) Dp , 

which implies 



12 



E 
and 

E [B[T\] 



E lH ( b )< L 

ben[T\ 



< n s {p+l)WTnL 2 Dp, 
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E la B <L 
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B[T] 

E X H(B)<L 
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+ E 


y~] ^a B >L 
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+ E 


2_j ^ a B>L 
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< n s (p + \)WTnL 2 Dp + E 



B[T) 

E la B >L 
5 = 1 



(23) 

Next, applying the Cauchy-Schwarz inequality to inequality 
(f79] >. we can obtain that 



f B[T] 

E aB 

,B=1 



which implies that 
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> LE 
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> L(E[B[T]]~n s (p+l)WTirL 2 Dp) 

where the first inequality follows from the Jensen's inequality 
and the last inequality follows from inequality ( l23l 

Sinc e the inequality holds for any L > 0. By choosing 

L ^ _ / E[B[T]] - 



2WTTrn s (p+l)pD 

"32 



, we can obtain that 



-gWT(p+l)yfcD>E[B[T\]. 



After substituting into d23l . we have 

"32 



E[A[T}] < 5K\og(n s p)[ x / — WT(p+l)^D 



IQkWT 
A 5 " 



p(logp) log(n s _p). 



